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Abstract 



This paper studies the capacity of the multi-antenna or multiple-input multiple-output (MIMO) secrecy 
channels with multiple eavesdroppers having single/multiple antennas. It is known that the MIMO secrecy 
capacity is achievable with the optimal transmit covariance matrix that maximizes the minimum difference 
^ ' between the channel mutual information of the secrecy user and those of the eavesdroppers. The MIMO 

secrecy capacity computation can thus be formulated as a non-convex max-min problem, which cannot be 
solved efficiently by standard convex optimization techniques. To handle this difficulty, we explore a relationship 
between the MIMO secrecy channel and the recently developed MIMO cognitive radio (CR) channel, in which 
the multi-antenna secondary user transmits over the same spectrum simultaneously with multiple primary users 
having single/multiple antennas, subject to the received interference power constraints at the primary users, or the 
so-called "interference temperature (IT)" constraints. By constructing an auxiliary CR MIMO channel that has 
the same channel responses as the MIMO secrecy channel, we prove that the optimal transmit covariance matrix 



to achieve the secrecy capacity is the same as that to achieve the CR spectrum sharing capacity with properly 
selected IT constraints under certain conditions. Based on this relationship, several algorithms are proposed to 



solve the non-convex secrecy capacity computation problem by transforming it into a sequence of CR spectrum 
sharing capacity computation problems that are convex. For the case with single-antenna eavesdroppers, the 
proposed algorithms obtain the exact capacity of the MIMO secrecy channel, while for the case with multi- 
antenna eavesdroppers, the proposed algorithms obtain both upper and lower bounds on the MIMO secrecy 
capacity. 
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Index Terms 

Cognitive radio, convex optimization, CR spectrum sharing capacity, interference temperature, multi-antenna 
systems, secrecy capacity. 

I. Introduction 

In the 1970s, Wyner introduced a secrecy transmission model in his seminal work [1] on information- 
theoretic secrecy. In this model, the secrecy transmitter sends confidential messages to a legitimate 
receiver subject to the requirement that the messages cannot be decoded by an eavesdropper. The 
information-theoretic study of the secrecy transmission problem has been continued and extended to 
many other channel models, including broadcast channels (BCs) [2], [3], multiple access channels 
(MACs) [4], [5], and interference channels (ICs) [6], [7]. Very recently, the secrecy capacity of the 
multi-antenna/multiple-input multiple-output (MIMO) channel has been characterized by Khisti and 
Wornell [8], and Oggier and Hassibi [9]. In their work, the MEMO secrecy channel with a single 
eavesdropper having multiple antennas is transformed into a degraded MIMO-BC, whose capacity is an 
upper bound on the secrecy capacity. It was shown in [8], [9] that this capacity upper bound is indeed 
tight for the Gaussian noise case, i.e., the exact secrecy capacity. However, this computable secrecy 
capacity cannot be extended to the general case of multiple eavesdroppers. Moreover, Liu and Shammai 
[10] also established the MIMO secrecy capacity by using the channel enhancement technique [11]. 
However, no computable characterization of the secrecy capacity was provided in [10]. 

On the other hand, cognitive radio is considered as an efficient technology to dramatically improve 
spectrum utilization, thus having great potential to solve spectrum scarcity problem. In a spectrum- 
sharing CR system, the CR user or the so-called secondary user (SU) is allowed to simultaneously 
transmit with the licensed primary user (PU) over the same spectrum, provided that the SU to PU 
interference level is regulated below a predefined threshold, which is also called the "interference 
temperature (IT)" constraint. The capacity achieving transmission problems under the IT constraint for 
the secondary users have been studied in [12], [13], and [14] for the CR MIMO point-to-point channel, 
the CR MIMO-MAC, and the CR MIMO-BC, respectively. Since the IT constraint is a linear function 
of the transmit covariance matrix, the capacity characterization problem for the CR MIMO channel can 
be formulated as a convex optimization problem, and is thus solvable via the standard interior point 
method [12]. It is worth noting that the system models of the secrecy channel and the CR channel are 
fairly similar in the sense that the secrecy and SU transmitters need to regulate the resultant signal power 
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level at the eavesdropper and PU, respectively, so as to achieve the goals of confidential transmission 
and PU protection, respectively. 

In this paper, we study the capacity computation problem for the general case of the MIMO secrecy 
channel with multiple eavesdroppers having single/multiple antennas. Based on the results in [8], [9], 
the related MIMO secrecy capacity can be obtained via optimizing over the transmit covariance matrix 
of the secrecy user to maximize the minimum difference between the mutual information of the secrecy 
channel and those of the channels from the secrecy transmitter to different eavesdroppers. It can 
thus be shown that the resulting capacity computation problem is a non-convex max-min optimization 
problem, which cannot be solved efficiently with standard convex optimization techniques. To handle 
this difficulty, we consider an auxiliary CR MIMO channel with multiple PUs having single/multiple 
antennas and the same channel responses as those in the MIMO secrecy problem. We next establish 
a relationship between this auxiliary CR MIMO channel and the MIMO secrecy channel by proving 
that the optimal transmit covariance matrix for the secrecy channel is the same as that for the CR 
channel with properly selected IT constraints for the PUs under certain conditions. Based on such a 
relationship, we transform the non-convex MIMO secrecy capacity computation problem into a sequence 
of CR capacity computation problems, which are convex and thus can be efficiently solved. For the case 
of single-antenna eavesdroppers, the proposed algorithms obtain the exact capacity of the associated 
MIMO secrecy channel, while for the case of multi-antenna eavesdroppers, the proposed algorithms 
obtain both the upper and lower bounds on the MIMO secrecy capacity. 

The rest of this paper is organized as follows. Section [U presents the system models and problem 
formulations for the CR MIMO transmission and the secrecy MIMO transmission. Section ITU describes 
the main theoretical results of this paper on the relationship between the secrecy capacity and the CR 
spectrum sharing capacity. Section [IV] studies the case of single-antenna eavesdroppers, and develops 
several algorithms to compute the MIMO secrecy capacity. Section |V] extends the results to the case 
of multi-antenna eavesdroppers to obtain the upper and lower bounds on the MIMO secrecy capacity. 
Section [Vjj presents some numerical examples. Finally, Section [VTI1 concludes the paper. 

Notation: Uppercase boldface and lowercase boldface letters are used to denote matrices and vectors, 
respectively. (S) H , tr(5), and \S\ denote the conjugate transpose, the trace, and the determinant of a 
matrix S, respectively. 1Z K denotes the vector space of K x 1 real vectors, and 1Z denotes the field of 
real numbers, i" denotes an identity matrix. E[-] denotes statistical expectation. | ■ | denotes the absolute 
value of a complex number. 
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II. System Model and Problem Formulation 

In this section, we present system models and problem formulations for the CR MIMO transmission 
and the secrecy MIMO transmission in the following two subsections, respectively. 

A. CR MIMO Transmission 

As shown in Fig. EJa), we consider a CR MIMO channel, where the SU transmitter (SU-Tx) is 
equipped with N transmit antennas, and the SU receiver (SU-Rx) is equipped with M receive antennas. 
The SU-Tx to SU-Rx channel is denoted by a N x M matrix H s . Moreover, there are K single-antenna 
PU receivers denoted by PUj, i = 1, ■ ■ • , K, and the channel from SU-Tx to PU, is denoted by the 
N x 1 vector hi. The received signal y at SU-Rx is expressed as 

y = Hfx + z (1) 

where x is the transmit signal vector at SU-Tx, and z denotes the noise vector at SU-Rx. The entries 
of the noise vector are independent circularly symmetric complex Gaussian (CSCG) random variables 
of zero mean and covariance matrix I. Since the SU shares the same spectrum with the PUs, there are 
K IT constraints imposed to the SU transmission, expressed as E[|h.f a;| 2 ] < I\, i — 1, ■ • • , K, where 
Tj denotes the tolerable IT limit for PUj. 

Consider the CR MIMO transmission problem, in which we determine the optimal transmit covariance 
matrix for SU-Tx to maximize the data rate subject to the transmit power constraint and the IT constraints 
for the K PUs. Mathematically, this problem can be formulated as [12] 

(PA): max \og\I + HfSH s \ 

$ 

subject to: tr(5) < P 

where S = K[xx H ] denotes the transmit covariance matrix at SU-Tx, and P denotes the transmit power 
constraint. Note that S is a positive semi-definite matrix such that (PA) is a convex problem and can 
be solved efficiently by the standard interior point method [15]. 

B. Secrecy MIMO Transmission 

As shown in Fig. \V[b), we consider a MIMO secrecy channel, where the secrecy transmitter (SC-Tx) 
is equipped with N transmit antennas, and the secrecy receiver (SC-Rx) is equipped with M receive 
antennas. Moreover, there are K single-antenna eavesdroppers. In accordance with the earlier introduced 
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CR MIMO channel, the channel response from SC-Tx to SC-Rx is denoted by H s , and the channel 
response from SC-Tx to the ith eavesdropper (EAj) is denoted by hi, i = 1, • • • , K. According to the 
secrecy requirement, the transmitted message W from SC-Tx should not be decoded by any of the 
eavesdroppers, i.e., H(W\yi) > r, Vi, with j/j denoting the received signal at EA,, and r denoting the 
secrecy transmit rate. According to the results in [8], [9], the secrecy capacity can be obtained by solving 
the following optimization problem 



where S denotes the transmit covariance matrix of SC-Tx, similar to that of SU-Tx in the CR case, 
and of denotes the variance of the zero-mean CSCG noise at EAj. 

We see that (PB) is a non-convex optimization problem since its objective function is the difference 
between two concave functions of S and thus not necessarily concave. Existing methods in the literature 
[8], [9], [16], [17] for the MIMO secrecy capacity computation is only applicable to the case of a single 
eavesdropper. However, these methods cannot solve the case with multiple eavesdroppers (PB) even 
for the case where each eavesdropper has a single antenna^. 

Remark 1: According to Fig. [Q it is easy to observe that the system models of the CR transmission 
and the secrecy transmission bear the similarity that they both need to control the received signal power 
levels at both PUs and eavesdroppers. However, note that (PA) guarantees that the interference power 
at each PU receiver is below the required threshold without considering the PU noise power, while for 
(PB), through the second term in the objective function, the confidential level at each eavesdropper 
is not only related to the received signal power from SC-Tx, but also related to the noise power at 
eavesdroppers. Therefore, one immediate question is whether there exists a relationship between these 
two systems such that we can solve the non-convex problem (PB) by transforming it into some form 
of (PA) that is convex and thus efficiently solvable. With this motivation, we first study the relationship 
between these two problems, and then propose corresponding algorithms to solve (PB). 

III. Relationship Between Secrecy Capacity and CR Spectrum Sharing Capacity 

In this section, we present main theoretical results of the paper on the relationship between the secrecy 
capacity computation problem (PB) and the CR spectrum sharing capacity computation problem (PA). 

'Problem (PB) in the case of multi-antenna eavesdroppers will be studied later in Section Ivl 



(PB) : max min log 1 1 + Hf SH S \ - log 1 + 




subject to: tr(5) < P 



DRAFT 



6 



While the developed relationship applies to both single-antenna and multi-antenna CR/secrecy channels, 
we are particularly interested in the multi-antenna case since it provides a general guidance for solving 
(PB). 

Proposition 1: For a given (PB), there exists a set of IT constraint values, Ti, i = 1, ■ ■ ■ , K, such 
that the resulting (PA) has the same solution as that of (PB). 

Proof: Please refer to Appendix [A] ■ 

Proposition Q] establishes the relationship between (PA) and (PB). To further investigate this rela- 
tionship, we define an auxiliary function of TjS as 

g(T u ,T K ):= max \I + H^SH S \ 

s 

subject to: tr(5) < P ( 2 ) 

h?Sh i <T u i = l,--- ,K. 
Note that the only difference between Problem © and (PA) is that the objective function in Problem 
© does not involve a logarithmic function of matrix determinant while that in (PA) does. As a result, 
Problem © is non-convex since its objective function is not concave in S. Also note that Problem © 
is equivalent to (PA) since they have the same optimal solution for S. Therefore, although Problem 
© is non-convex, its optimal solution can be obtained via solving the convex counterpart (PA). 
Proposition 2: (PB) is equivalent to the following optimization problem: 

t?(t -r \ 9(^1, • • • ,r K ) 

max mm Fi[T u ■■■ ,T K ) := . (3) 

"lj— Fk * 1 + 1 i/ cr i 

Proof: Please refer to Appendix [Bj ■ 
Proposition [2] establishes the relationship between (PB) and the auxiliary function g(r l7 ■ • • , IV) that 
is related to (PA). The equivalence between Problem © and (PB) means that by solving the optimal 
T/S in Problem ©, we could solve an optimal S given g(Ti, ■ ■ ■ ,T K ) is an embedded optimization 
problem over S inside Problem ©. Such an optimal S is also the solution for (PB), for which the 
explanation is given in Appendix [Bj 

Problem © can be solved by utilizing an important property of g(T 1 , ■ ■ ■ , IV) described as follows: 

Proposition 3: The function g(Ti, ■ ■ ■ , Tk) is a concave function with respect to T±, ■ ■ ■ , IV, and 

7*(ri, • • • , rV) := ^^^^ = + Hf SUH S \, i = 1, • - - , K (4) 

where and /4 are the optimal solution of (PA) and the corresponding Lagrange multiplier (the 
dual solution) with respect to the ith IT constraint, respectively. 
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Proof: Please refer to Appendix O ■ 
Note that from Proposition [3l it follows that the gradient of g(Ti, • • ■ , IV) in © can be obtained by 
solving (PA) via the Lagrange duality method, which completes the equivalence between (PA) and 
(PB) via the intermediate problem ©. At last, we have 
Proposition 4: Problem © is a quasi-concave maximization problem. 

Proof: Please refer to Appendix [D] ■ 
Proposition H] suggests that Problem © can be solved by utilizing convex optimization techniques, for 
which the details are given in the next section. 

IV. Algorithms 

In this section, we present the algorithms to compute the MIMO secrecy capacity by exploiting the 
relationship between the secrecy transmission and the CR transmission developed in Section [TTTJ The 
algorithm for the general case of the MEMO secrecy channel with multiple eavesdroppers is presented 
first. Two reduced-complexity algorithms are next presented, one for the special case with one single 
eavesdropper, and the other for the special case with a single-antenna secrecy receiver, i.e., the multiple- 
input single-output (MISO) secrecy channel. 

A. General Case 

In this subsection, we present the algorithm for (PB) in the general case of MIMO secrecy channels 
with multiple eavesdroppers. According to Propositions [2] and HI (PB) is equivalent to the quasi-concave 
maximization problem ©. Thus, we instead study Problem © since it is easier to handle than (PB). 

According to [15], a quasi-concave maximization problem can be reduced to solving a sequence of 
convex feasibility problems. Thus, Problem © can be further transformed as 

max t 

subject to : g(T u ■■■ , T K ) > t(l + iy<7 2 ), i = 1, • • • , K. 
Let t* be the optimal solution of Problem ©. Clearly, t* is also the optimal value of Problem ©. If 
the feasibility problem 



max 

rv-.rv (6) 
subject to : g(T u ■■■ ,T K )>t(l + Ti/a?), % = 1, • • • , K 
for a given t is feasible, then it follows that t* > t. Conversely, if Problem (J6]) is infeasible, then t* < t. 
Therefore, by assuming an interval [ 0, t ] known to contain the optimal t*, the optimal solution of 
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Problem © can be found easily via a bisection search. Note that a suitable value for t can be chosen 
as (7(00, • • • , 00) from ©. 

We next solve the feasibility problem © by a similar method discussed in [18]. It is worth noting 
that the feasibility problem © can be viewed as an optimization problem. The Lagrangian of Problem 
© can be written as 

K 

L ({ui}, r 1; • ■ ■ , r K ) = (</(ri, ■■■,r K )-t(i + ry<r?)) (7) 

i=l 

where is the non-negative dual variable for the ith constraint, and {z^} denotes v\, • ■ ■ ,vk- The 
corresponding dual function is then defined as 

K 

/ (W)= max V)z/ i (y(r 1 ,...,r ii .)-t(i + r i / < 7?)). (8) 

2=1 

Due to its convexity, Problem © can be transformed into its equivalent dual problem as 

min / (M) (9) 

{^} 

and the duality gap between the optimal values of Problem © and Problem © is zero if Problem © 
is feasible. 

Since it is known from Proposition [3] that function g(T 1 , ■ ■ ■ , T K ) is concave with respect to {r 1; • • • , 
IV}, Problem ® can be solved via a gradient-based algorithm. According to Proposition [3] the gradient 
of function g(T 1 , ■ ■ ■ ,T K ) can be obtained by solving (PA). Furthermore, since function fo({vi}) is 
convex with respect to {^}, Problem © can be solved by a subgradient-based algorithm, such as the 
ellipsoid method [15]. Similar to Lemma 3.5 in [18], Problem © is infeasible if and only if there exist 
{ui} such that fo({vi}) < 0. Using this fact along with the subgradient-based search over {^}, the 
feasibility problem © can be solved. In summary, the algorithm for Problem © with a target accuracy 
parameter e is listed as follows: 
Algorithm 1: 

. Initialization: t min = 0, t max = i. 

• Repeat 

- t <- |(t min + t max ). 

- Solve the feasibility problem ©. If Problem © is feasible, t mm t; otherwise, t max <— t. 

- Stop, when t max - t min < e. 

• The optimal value of Problem © is taken as t mm . 
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B. Single-Eavesdropper Case 

We now consider a special case of (PB), where there is only one single eavesdropper in the secrecy 
channel, and propose a simplified algorithm over Algorithm 1 for the corresponding (PB). 

Consider first the counterpart CR transmission problem (PA). For the single-PU case, (PA) can be 
rewritten as 

(PC): max \og\I + HfSH s \ 

s 

subject to: tr(S) < P 

h H sh < r 

where h denotes the channel from SU-Tx to the single PU, and V is the corresponding IT limit for the 
PU. On the other hand, for the single-eavesdropper case, the secrecy transmission problem (PB) can 
be rewritten as 

/ h H Sh\ 

(PD): max log|I + JEff SJEf a | -log(l + =— ] 

S \ a z / 

subject to: tr(5) < P 

where h denotes the channel from SC-Tx to the single eavesdropper, and a 2 denotes the variance of 
the noise at the eavesdropper. Following Proposition [2l (PD) is equivalent to the optimization problem 

max F(r) := (10) 

where the function g(Y) is the single-PU counterpart of that in ©. Moreover, it is evident from 
Proposition H] that function F(T) is quasi-concave and the optimization problem (fTOl ) is a quasi-concave 
maximization problem. 

Lemma 1: The sufficient and necessary condition for T* to be the optimal solution of Problem (fTOl ) 

is 

A 1 



1 (r)(l + T*/a 2 ) = -g(T*) (11) 
a 2 - 



dg{T) 



where 7(1") :-- 9T . 

Proof: Please refer to Appendix [0 ■ 
Based on Lemma [Q (PD) can be solved via the equivalent problem (flOl) by the bisection method 

summarized as follows: 

Algorithm 2: 

. Initialization: T min = 0, T max = f . 
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. Repeat 

- r<- |(r min + r max ). 

- Solve Problem © for the single-PU case, and compute 7(T). If 7(r)(l + Y/a 2 ) > ^sg(T), 
r min ^_ r . ot herwise, T max <- r. 

- Stop, when T max - T min < e. 

• The optimal solution of (PD) equals that of (PC) with the converged T. 
Note that in the above algorithm, f = h H S Q h and S is the optimal solution of (PC) without the 
interference power constrain^]. 

Algorithm 2 searches the optimal T according to its gradient direction, and thus avoids solving the 
sequence of feasibility problems in Algorithm 1. Therefore, Algorithm 2 is much simpler than Algorithm 
1. However, since the general case of (PB) has multiple variables I^s, this gradient-based algorithm 
cannot be applied to the general case. 

Remark 2: Similar to Proposition [TJ a dual relationship between the secrecy transmission (PD) and 
the CR transmission (PC) in the case of a single eavesdropper/PU can be described as follows. For a 
given (PC), there is a parameter a, such that (PD) with the noise variance a 2 at the eavesdropper has 
the same solution as that of (PC). This property can be proved by combining Lemma[Hand Proposition 
|2] This proof is thus omitted for brevity. 

C. Single-Antenna SC-Rx Case 

We now turn our attention to another special case of the secrecy channel where SC-Rx is equipped 
with a single receive antenna, i.e., the MISO secrecy channel. Same as the MIMO secrecy case, each 
eavesdropper is still assumed to have a single antenna. For notational convenience, (PA), (PB), (PC), 
and (PD) in the case of single- antenna SC-Rx are denoted correspondingly by PA-s, PB-s, PC-s, and 
PD-s. 

The problem PC-s has been studied in [12]. In [12], it was shown that the optimal transmit covariance 
matrix for the CR MISO channel is a rank-one matrix, and a closed-form solution for the optimal transmit 
beamforming was presented. The problem PD-s has been studied in [6], [16], where it was shown that 
the optimal transmit covariance matrix for the secrecy MISO channel is also a rank-one matrix, and 
based on the generalized eigenvalue decomposition, a closed-form solution for the optimal transmit 
beamforming was provided. 

2 When r > f , the value of g(T) is constant regardless of V. Thus, the optimal solution of Problem ilOi satisfies T* < f . 
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Consider PA-s, in which there are multiple PUs each having a single receive antenna. To the authors' 
best knowledge, no closed-form solution exists for such a case. Nevertheless, due to its convexity, this 
problem can be solved via a standard interior point algorithm. By using a similar method to that in 
[12], it can be shown that the optimal transmit covariance matrix for PC-s is also a rank-one matrix. In 
contrast, for PB-s, due to its non-convexity, there is no existing method in the literature to solve this 
problem. However, since PB-s is a special case of (PB), we can apply Algorithm 1 to efficiently solve 
this problem . 

Next, by exploiting the special structure of PB-s, we provide a simplified algorithm over Algorithm 
1. First, we rewrite PB-s as 



(PB-s) : maxmin Fi(S) :- 



1 + h* Sh s 



S ' I • (h?Stn) aj (12) 

subject to: tr(S) < P 

where the JV x 1 vector h s denotes the channel from SC-Tx to the single antenna SC-Rx. Unlike the 
general case of (PB) where only its transformed problem in © is a quasi-concave problem with respect 
to TjS, PB-s itself is a quasi-concave problem with respect to S due to the following proposition. 
Proposition 5: F^S) is a quasi-concave function for i = 1, . . . , K. 

Proof: Please refer to Appendix |FJ ■ 
Thus, PB-s can be transformed into the following equivalent problem 

max t 
S,t 

subject to: tr(5) < P (13) 

l + hf5h.>t(l + ^^),i = l,-.. ,K 

where t is a positive variable. For the fixed t, all the constraints in the above problem are linear matrix 
inequalities over S, and thus the corresponding feasibility problem (similarly defined as ©) can be 
viewed as a semi-definite programming (SDP) feasibility problem. Correspondingly, the optimal value 
of t can be obtained by a bisection search algorithm. 

Compared with Algorithm 1, the algorithm for Problem (fl"3T) is much simpler, since the SDP feasibility 
problem can be solved via high-efficiency interior point methods, while the feasibility problem © in 
Algorithm 1 can only be solved through a general gradient-based algorithm. Moreover, according to 
Proposition [Q we can find a set of parameters I^s such that the corresponding PA-s has the same 
solution of PB-s. Since the optimal solution of PA-s is known to be a rank-one matrix [12], so is the 
optimal solution for PB-s. 
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V. Multi-antenna Eavesdropper Receiver 

In this section, we extend our results to the case with multi-antenna eavesdroppers. We assume 
that each eavesdropper is equipped with iV e receive antennas, and the channel from SC-Tx to the ith 
eavesdropper receiver is denoted by Hi of size N x N e . Similar to (PB), the MIMO secrecy capacity 
in the multi-antenna eavesdropper case can be obtained from the following optimization problem [9] 



^PE) : max min log | / + Hf SH S | - log | / + Hf SHi | (14) 

S i 



subject to: tr(S) < P. (15) 

To the best knowledge of the authors, there is no existing solution in the literature for the above 
problem. In the following, we derive the upper and lower bounds on the MIMO secrecy capacity in 
the multi-antenna eavesdropper case based on the relationship between the secrecy transmission and the 
CR transmission. 

A. Capacity Lower Bound 

First, we have the following lemma: 

Lemma 2: If for any i, i G {1, • • • .K}, ix(HfSH t ) < T u we have \I + HfSH,] < (1 + ^) L , 
where L = min(A^ e , N). 

Proof: Please refer to Appendix [Gl ■ 
Similar to Proposition [2l from Lemma [2l the following proposition holds: 

Proposition 6: The optimal value of (PE) is lower-bounded by that of the following optimization 
problem 

p f-r -r \ SK^I) • • • j Pfc) n „ 
max mm F i (T 1 , ■ • ■ , T K ) := — — f — (16) 



ri,-,r 



K 



where the function g(Ti, ■ ■ ■ ,T K ) is defined as 

g(T u ■ ■ ■ , T K ) := max|/ + H^SH S \ 

s 

subject to: tr(S) < P ( 17 ) 

\x{Hf SH^ < r h i = l,---,K. 
Problem (TT6l) can be solved by the gradient-based method similar to Algorithm 1. Accordingly, the 

lower bound on the MIMO secrecy capacity is obtained. Note that this capacity lower bound is tight 

when N P = 1 and thus L = 1. 
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B. Capacity Upper Bound 

In the multi-antenna eavesdropper case, the signals received at different antennas of each eavesdropper 
are jointly processed to decode the contained secrecy message. Therefore, a straightforward upper bound 
on the secrecy capacity in this case is obtained by assuming that the signals at different antennas 
of each eavesdropper are decoded independently. Suppose that h^j is the jth column of the matrix 
Hi,j = !,■■•, N e , then the upper bound on the secrecy capacity can be obtained as 



The above problem is the same as (PB) with the number of single-antenna eavesdroppers equal to 
N e K, and thus can be solved by Algorithm 1. 



In this section, we provide several numerical examples to illustrate the effectiveness of the proposed 
algorithms in computing the secrecy channel capacity under different system settings. For all examples, 
we consider a MIMO secrecy channel with M = N = 4. The elements of the matrix H s and the vectors 
hiS (or the matrices His, in the multi-antenna eavesdropper case) are assumed to be independent CSCG 
random variables of zero mean and unit variance. Moreover, the noise power at each eavesdropper 
antenna is chosen to be one, and the transmit power of the secrecy transmitter, P, is defined in dB 
relative to the noise power. 

A. MIMO Secrecy Capacity with Two Single-Antenna Eavesdroppers 

In this example, we consider a MIMO secrecy channel with K = 2 single- antenna eavesdroppers. 
Fig. [2] plots the secrecy capacity of this channel obtained by Algorithm 1, where the transmit power 
ranges from dB to 10 dB. Moreover, a reference achievable secrecy rate of this channel is obtained by 
the Projected-Channel SVD (P-SVD) algorithm in [12]. In this algorithm, the channel H s is projected 
into a space, which is orthogonal to hi and h 2 , and thus the secrecy signals cannot be received by the 
eavesdroppers. It is easy to observe from Fig. [2] that the secrecy rate obtained by P-SVD is less than the 
secrecy capacity obtained by Algorithm 1 . Moreover, from Proposition HI it is known that the function 
Fi(Ti, T 2 ) is a quasi-concave function, and thus the function minj =12 F^Ti, L 2 ) is also a quasi-concave 
function. In Fig. |3] we plot the value of this function for P = 5 dB. It is observed that this function is 
indeed quasi-concave. 



max min log \ I + H Q SHJ — log ( 1 + 

S V 




(18) 



subject to: tr(S) < P. 



VI. Numerical Examples 
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B. MIMO Secrecy Capacity with One Single-Antenna Eavesdropper 

In this example, we apply Algorithm 2 to compute the secrecy capacity of a MIMO channel with 
one single-antenna eavesdropper. As shown in Fig. HI the secrecy capacity obtained by Algorithm 2 is 
larger than the achievable secrecy rate obtained by the P-SVD algorithm. Moreover, it is verified that 
function F(T) defined in (flOl) is indeed quasi-concave in Fig. [5] for P = 5 dB. 

C. MIMO Secrecy Capacity with One Multi-antenna Eavesdropper 

In this example, by applying the methods discussed in Section |Vj we show in Fig. |6] the lower and 
upper bounds on the MIMO channel secrecy capacity with a single eavesdropper using N e = 2 receive 
antennas. From the capacity lower bound, we obtain a feasible transmit covariance matrix and thus 
a corresponding achievable secrecy rate, shown in Fig. [6] and named as "Achievable Secrecy Rate". 
Moreover, the achievable secrecy rate by the P-SVD algorithm is also shown for comparison. 

VII. Conclusion 

In this paper, we have disclosed the relationship between the multi-antenna CR transmission problem 
and the multi-antenna secrecy transmission problem. By exploiting this relationship, we have transformed 
the non-convex secrecy capacity computation problem into a quasi-convex optimization problem, and 
developed various algorithms to obtain the optimal solution for different cases of secrecy channels. 
Although the proposed method cannot obtain the exact secrecy capacity for the more complicated 
multi-antenna eavesdropper case, it can be applied to compute the upper and lower capacity bounds. 

Appendix 

A. Proof of Proposition^ Proposition Q] can be proved by contradiction. For the fixed er^s, suppose 
that the optimal solution of (PB) is S D . Define f j = hf S D hi, i — 1, . . . , K. If the optimal solution of 
(PA) with Yi = T it Vi, denoted by S Q , satisfies log 1 1 + Hf S H S \ > log 1 1 + Hf S H S \, then S D is 
a better solution for (PB) than S , which contradicts the preassumption that S is the optimal solution 
of (PB). Then there must be log \I + H^S H S \ < log \I + HfS Q H s \, which means that S Q is also 
the optimal solution of (PA), with Ti = hf S hi, i — 1, . . . , K. Proposition Q] thus follows. 

B. Proof of Proposition^ It is easy to observe that (PB) can be re-expressed as 

. \I + HfSH s \ 

max mm 77 

S i 1 + hf Shi/a? (19) 

subject to: tr(5) < P. 



DRAFT 



15 



Suppose that S a is the optimal solution of Problem ( fT9l ) and (PB). Define T a :— \ I + Hf S Q H S \ and 
f i := hfS hi,i = 1, • • • ,K, then the optimal objective value of Problem (fT9l) is F = min^T D /(l + 

r i),---, T /(i+r K )y 

Suppose that the optimal solution S a of Problem © with = fj, Vz, satisfies \I + S H S \ > T a , 
then S is a better solution for Problem (fl9l) than So, which contradicts the preassumption that S is 
the optimal solution of Problem (fT9l) . Therefore, we have T Q = (/(Fi, • • • , T K ). Thus, F is achievable 
for Problem © with the particular choice of r, = f is Vi. 

Furthermore, suppose that r\s are the optimal solutions of Problem ©, and the corresponding optimal 
objective value is F. For Problem © with Tj = Tj, suppose that the optimal solution is S. We can 
prove that F < F by contradiction: If F > F, S is a better solution for Problem (fT9l) than ,S , which 
contradicts the preassumption that So is the optimal solution of Problem (fT9l) . As such, we see that 
F is not only achievable for Problem ©, but also the optimal value of Problem © with the optimal 
solutions given as S = S Q and fj = hf S Q hi, Vi (Note that S is a hidden design variable for Problem 
©.)■ Proposition [2] thus follows. 

C. Proof of Proposition \3\- We first study several important properties of Problem © that is known 
to be an equivalent problem of (PA). Considering (PA) first, its Lagrangian function can be written 

as 

K 

L X {S, A, {m}) = log \I + HfSH s \ - A(tr(5) - P) - ^ fi^hfSh, - I\) (20) 

i=l 

where A and /ij are the non-negative Lagrange multipliers/dual variables with respect to the transmit 
power constraint and the interference power constraint at PUj, respectively. Since (PA) is a convex 
optimization problem, the Karush-Kuhn-Tucker (KKT) conditions [15] are both sufficient and necessary 
for a solution to be optimal, and solving (PA) is equivalent to solving its dual problem 

min maxL 1 (S', A, {yUj}). (21) 

On the other hand, the auxiliary problem © is non-convex due to the fact that its objective function 
is not concave. In general, the KKT conditions may not be sufficient for a feasible solution to be optimal 
when we have a non-convex optimization problem. However, we prove in the following lemma that this 
is not the case for Problem ©. 

Lemma 3: With Problem ©, the KKT conditions are both sufficient and necessary for a solution to 
be optimal. 
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Proof: The necessary part of Lemma [3] is obvious even for a non-convex optimization problem 
[15]. The sufficient part of Lemma [3] can be proved via contradiction as follows. The Lagrangian of 
Problem ([2]) can be written as 

K 

L 2 (S, 5, { lt }) = \I + HfSH s \ - 5(tr(S) - P) - ]T ^(hfShi - T z ) (22) 

i=X 

where 5 and 7, are the non-negative dual variables with respect to the transmit power constraint and 
the interference power constraint at PUj, respectively. We first list the KKT conditions of Problem © 
as follows: 

K 

\I + HfSH s \H s (I + HfSH^Hf = SI + J2 H h i h f (23) 

i=l 

5(tr(S) - P) = (24) 
^(hfShi-Ti) = 0, 1 = l,--- ,K. (25) 

Suppose that S^°\ 5^°\ and 7-°' ) are a set of primal and dual variables that satisfy the above KKT 
conditions, and the corresponding optimal value of Problem (O is C^°\ 
The KKT conditions of (PA) are expressed as 

K 

H S (I + HfSH^Hf = \I + Y J frhih? (26) 

i=l 

A(tr(5) - P) = (27) 
fi i (hfSh l -T l ) = 0, i = l,... (28) 

Suppose that S^ 1 -*, and /i^ are the optimal primal and dual variables that satisfy the above KKT 
conditions, and the corresponding optimal value of (PA) is C^ 1 ' . Note that since (PA) is convex, the 
KKT conditions are both necessary and sufficient. 

If (|23"l)-(l25]l are not sufficient such that log(C (0) ) ^ i.e., 5 (0) ^ we can choose 

S = S (0) (29) 
\ = 5 {0) /\I + HfS m H s \ (30) 
m = ^ 0) /\I + HfS^H t \, i = l,-..,K (31) 

for (PA), which clearly also satisfy the KKT conditions of (PA). Given the sufficiency of the KKT 
conditions for (PA), 5 (0) is also optimal for (PA) based on ([29]) such that log(C (0) ) = C (1) , which 
contradicts our assumption that log(C^ ^) 7^ Lemma [3] thus follows. ■ 
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Essentially, it is due to the equivalence between the non-convex Problem © and the convex (PA) 
that Lemma [3] holds. From Lemma it follows that the duality gap between Problem © and its dual 
problem, defined as 

D = min max L 2 (S, 5, {7*}), (32) 

<5,{7i} S 

is zero, i.e., g(T 1 , ■ ■ ■ , T K ) = min^^} maxg L 2 (S, 5, {%})■ As such, from (|22~1) we have 

dg(T 1 ,.~ t r K ) = dD = (Q) 

ar, ar^ 7l ' ' ' p ; 

Combining (1311) and (1331) . the latter part of Proposition [3] thus follows. 

Now we prove the concavity of g(Ti, ■ ■ ■ , T K ). For the function g(q), where q := [r 1; ■ ■ • , T K ] T E 
1Z+, its concavity can be verified by considering an arbitrary line given by q = x + tv, where x E 71+ , 
v E 1Z K , t E 1Z + , and x + tv E 1Z+ [15]. In the sequel, we just need to prove that the function g(x + tv) 
is concave with respect to t. Moreover, if the ith IT constraint is not active for Problem ©, we have 
7i = from the KKT condition such that the concavity holds. To exclude the above trivial case, we 
assume that all K IT constraints are active for Problem © in the following. 

Define 

f 2(8,1/1, ■ ■ ■ ,1k) '= max L 2 (S, 5, 71, • • • ,i K ) (34) 

iJ 

as the dual function of Problem ©. Let s be the subgradient of f 2 (S, 71, ■ • • ,1k)- According to the 
definition of subgradient, the subgradient at the point [5, 71, ■ • • ,1k] satisfies 

f2(S, 71, • • • , Ik) > /2 ft li, • • • , Ik) + ([£, li, • • • , 7k] - [5, 7i, ' • • ,1k])- s , (35) 

where [5, 71, • • ■ , ik\ is another arbitrary feasible point. 

Lemma 4: The subgradient s of function f 2 (5, 71, • ■ ■ , 7k) at point [5, 71, • • ■ , 7V] is [P — tr(5), i\ — 
h 1 ^ Shi, ■ ■ ■ , — h^STiR-], where S is the optimal solution of Problem (1341) at this point. 

Proof: Let S 1 be the optimal solution of Problem (1341) with 5 = 5 and 74 = 7^ i = 1, • • • , i^. Thus, 
we have 

K 



/ 2 (5, 71, • • • , Ik) = f - 5(tr(5) - P) - ^ 7i (/if Sfc - I\) (36) 

i=l 
if 

> f - 5(tr(5) -P)-J2 li( h f~ Sh i ~ T d ( 37 ) 



i=l 
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K K 



= r- ~5(tv(S) - P) - ^jiihfShi - Ti) + %(S) - P) + J>i(hf Sh { - T t ) 

i=i i=i 

K 

- S(tr(S) - P) -Y,li{h?Shi - Ti) (38) 

i=i 

K 

= / 2 (*,7i, ■■■n K ) + (tr(S) - P)(8 -8) + ^(hfShi - r,)(7 4 - 7*) (39) 

i=i 

where f = |J + HfSH s \ and f = |J + HfSH s \. According to (1391 , we have Lemma @J ■ 
According to Lemma [3l Problem © is equivalent to its dual problem 

min / 2 (5,7i, • • • ,7k) (40) 

o,7i >■■■ ,7-ff 

where / 2 (5, 71, • • • ,7k) is convex. We next consider Problem © with parameters P, Ti, • • • , IV, 
denoted as Problem I. Assume that <^,7i , • • • ,7^ are its primal and dual optimal solutions. 
Moreover, we have another form of Problem © with parameters P, Ti + tvi, ■ ■ • , Tk + tvx, denoted as 
Problem II, where t is a positive constant and Vi is a real constant. Assume that S^ 2 \ 8^ 2 \ 7^, ■ • • , 7^ 
are the primal and dual optimal solutions of Problem II. According to (|34|) . we can write the dual 
function of Problem II as 

K 

/ 2 n (<5,7i, • • • ,1k) := max |J + iff 5/f s | - 5(tr(S) - P) - ^7i( h f Sh i " r * " tv i) ( 41 ) 

To solve Problem II, we apply the subgradient-based algorithm to search the minimum of its dual func- 
tion /^(S, 7i> " ' > Ik) along the subgradient direction. Suppose that we start from the point [5^\ 7^ , • ■ ■ , 
7^]. Based on Lemma H one valid subgradient of f 2 (5 } 71, ■ ■ • ,7k) at this point is 

[O^i + tUi-hfSWfci,-- - ^jf + tojf-fef^ 1 ^] = [0,t«i,--- (42) 

where (@2]) is due to the KKT condition of Problem I: rf ) -hf S {1) hi = given 7^ > 0, Vi. Moreover, 
according to (1351) . we have 

/?(^\7?\"^7^ (43) 

where is the subgradient at the point [5^\ ^\ ■ ■ ■ ,7^]. Since 6^, ^[ 2 \ ■ ■ ■ ,7^ are the dual 
optimal solutions of Problem II, we have /"O^jTi , • ■ • ,7k ) — /^O^jTi? ' ' ' ?7k' ) )- Combining 
this with (O and (|43l . we have 

k k 

J>f^i< J>fV (44) 

8=1 1=1 
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Thus, 



K K 



11 'Vi 



i=l i=l 



< ^2li 1] Vi, given t > 0. (45) 



Moreover, according to Lemma [3] and (|22|) . we have 

*<5+^ = £>. (46 ) 

i=l 

Note that 73 is the Lagrange multiplier of Problem © with respect to the ith IT constraint. With a 
different IT threshold, i.e., a different value of t, jiS are not necessarily the same, and thus 7iS can be 
viewed as implicit functions of t. Combining (l45l) with (l46l) . it is easy to observe d9( - x + tv ^ decreases 
with the increase of t since the derivative change over t is given as Ylf=\ n v i ~~ Si=i li^ v i — 0> i- e -> 
the second order derivative of function g(x + tv) over t is negative on an arbitrary line x + tv in the 
feasible region. Therefore, g(q) is concave. Proposition [3] thus follows. 

D. Proof of Proposition^ The quasi-concavity is define as follows [15]: 
Definition 1: A function / : 1Z K — > 72. is called quasi-concave if all its sublevel sets 

S a = {xe domf\f(x) > a} (47) 

for a E 71, are convex sets. 

According to Proposition [3] gfix, ■ ■ ■ , IV) is a concave function of TjS. Therefore, the a-sublevel set 
of Fi(Ti,--- ,T K ) 

s a = \q 9{T ^' r ; T 2 K) > <*} = {qWi, ■■■,r K )> a (i + r 4 /of )} (48) 

l 1 + 1 i/a/ J 

is a convex set for any a, and thus the function Fi(Tx, ■ ■ ■ ,Tk) is a quasi-concave function. Since 
the objective function of Problem © is the minimum of K quasi-concave functions, Fi(Ti, ■ ■ ■ , IV), 
% = 1, • • • , i^, it is still quasi-concave [15]. Proposition |4] thus follows. 

E. Proof of Lemma [7} The optimality condition of Problem (flOl) is 

ggg) _ 7(r)(i + r/a 2 )-^(r) _ 

(i + r/^) 2 1 } 

Since the above optimality condition is a necessary condition for any unconstrained smooth optimization 
problems regardless of its convexity [15], the necessary part of Lemma \T\ follows. 

We next prove the sufficient part of this lemma by contradiction. We first present a property of 7(T) 
as follows. 

Lemma 5: j(T) is a non-increasing function for T > 0. 
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The proof of Lemma [5] is similar to that of Proposition @] and thus is omitted here. Suppose that 
there are two solutions and T^, both of which satisfy the condition in (fTTj) . Furthermore, without 
loss of generality, we assume rW < T^. Therefore, we have 

g(r®) > g(T®). (50) 



According to (1331) . we have 

r(2)_rw - 7l y ' 

Thus, 

0(r( 2 )) - r^ 7 (r( 2 )) > <?(r«) - r« 7 (r( 2 )) > <?(r«) - r« 7 (r«) (52) 

where the second inequality is due to the fact that 7 (r( 1 )) > 7 (r( 2 )) from Lemma [51 
Since both solutions satisfy the necessary condition (fTTT) . we have 

7(r (1) ) = ^(^(r«)- 7 (rW)r«) (53) 
7(r (2) ) = ^(^(r( 2 ))- 7 (r( 2 ))r( 2 )). (54) 

From (1521) and 7 (r^^) > 7 (r^ 2 ^), it is easy to observe that (1531) and (l54l) hold simultaneously if and 
only if 7 (r^) = 7 (r^ 2 ^) = 7. Thus, we have 
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(T^)=g(T^)- 1 (T^ -T^). (55) 



Since < T^ 2 ) and g(T^) > g(r( 2 )), we further derive 7 < 0, which contradicts the fact that the 
Lagrange multiplier 7 > 0. As such the solution of ([111) is unique, which implies the sufficiency given 
the already proven necessity part. Lemma Q] thus follows. 



F. Proof of Proposition [3} Similar to the proof given in Appendix ITJ1 the a-sublevel set of F^S) 

1 + hfSh s 
1 + (hfSh^/ai 

{S\l + fcf Sh s > a(l + {hfShi)/al)}. (57) 



S a = {S\ 7 i r^>a} (56) 



is a convex set. Thus, Fi(S) is a quasi-concave function. 
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G. Proof of Lemma |3t We have 

\I + HfSHi\ = \I + UfAiUi\ = \I + Ai| 

where Hf SHi := Uf A;Z7j is the eigenvalue decomposition. Since \x(Hf SHi) 
tr(HfSHi) < Ti it follows that 

tr(Ai) < IV (59) 

Combining (|58l ) and (|59l) and denoting L = min(/V e , N), we have 

\I + H?SH t \ < \l+y-l\ = (1 + ^) L (60) 

where the inequality is obtained by solving the following problem: max tr ( A .)< r . |J+Aj|. Lemma [2] thus 
follows. 
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SU-Tx 




SC-Tx 



(b) 

Fig. 1. The comparison of two system models: (a) the CR MIMO channel with K single-antenna PUs; and (b) the secrecy 
MIMO channel with K single-antenna eavesdroppers. 
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Fig. 2. Comparison of the secrecy capacity by Algorithm 1 and the secrecy rate by the P-SVD algorithm in [12] for 
M = N = 4 and K = 2 single-antenna eavesdroppers. 
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Fig. 3. The value of the function min^i^ Fi(T\, I^) for M = N = 4, K = 2 single-antenna eavesdroppers, and P = 5 
dB. 
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4.5 




Transmit power (dB) 

Fig. 4. Comparison of the secrecy capacity by Algorithm 2 and the secrecy rate by the P-SVD algorithm in [12] for 
M = N = 4 and K = 1 single-antenna eavesdropper. 
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Fig. 5. The value of the function F(T) for M = N = 4, K = 1 single-antenna eavesdropper, and P = 5 dB. 
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-□- Secrecy Rate by P-SVD 
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. 6. Comparison of the lower and upper bounds on the secrecy capacity and two achievable secrecy rates for M = N = 4, 
= 1 eavesdropper with N e = 2 receive antennas. 
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